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ABSTRACT 

We examine the consistency of Hi properties with the unification scheme for radio 
galaxies and quasars, and any correlation with the symmetry parameters for a sample 
of CSS and GPS sources. In our sample, 15 out of 23 galaxies exhibit 21-cm absorption 
as against 1 out of 9 quasars, which is broadly consistent with the unification scheme. 
Also there is a tendency for the detection rate as well as the column density for galaxies 
to increase with core prominence, f c , a statistical indicator of the orientation of the 
jet axis to the line of sight. This can be understood in a scenario where radio sources 
are larger than the scale of the circumnuclear Hi disk so that the lines of sight to the 
lobes at very large inclinations do not intersect the disk. The sources in our sample 
also exhibit the known anti-correlation between Hi column density and source size. 
This suggests that small linear size, along with intermediate values of core prominence, 
is a good recipe for detecting 21-cm absorption in CSS and GPS objects. Some of the 
absorption could also be arising from clouds which may have interacted with the radio 
jet. The Hi column density and velocity shift of the primary absorption component, 
however, show no dependence on the degree of misalignment and the separation ratio 
of the radio sources. 

Key words: galaxies: active - galaxies: evolution - galaxies: nuclei - galaxies: ab- 
sorption lines - radio lines: galaxies - quasars: general 



1 INTRODUCTION 

An understanding of the distribution and kinematics of the 
different components of the circumnuclear gas in an active 
galactic nucleus (AGN) is important both for studying the 
anisotropy in the radiation field and thereby testing the uni- 
fication scheme (cf. Barthel 1989; Urry & Padovani 1995), 
and for understanding the fuelling of the radio activity. At 
radio wavelengths, the ionized component of this gas may 
be probed via radio polarization measurements of compact 
components in sources which reside within the dense inter- 
stellar environments of the parent galaxies. Such sources are 
the radio cores, and the compact steep-spectrum (CSS) and 
gigahertz peaked-spectrum (GPS) sources. The CSS sources 
are defined to be those with a projected linear size <15 
kpc (H =71 km s" 1 Mpc" 1 , fi m =0.27, fi A =0.73, Spergel 
et al. 2003) and having a steep high-frequency radio spec- 
trum (a>0.5, where S(v) oc u~ a ). The GPS sources are more 
compact and have sizes <1 kpc (O'Dea 1998). An important 
way of probing the atomic gas on sub-galactic scales is by 
studying 21-cm Hi absorption towards the compact compo- 
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nents of CSS and GPS objects, or the compact radio nuclei 
of the larger objects (e.g. van Gorkom et al. 1989; Conway 
& Blanco 1995; Peck et al. 2000; Pihlstrom 2001; Vermeulen 
et al. 2003; Pihlstrom et al. 2003; Gupta et al. 2006). 

The Hi absorption lines associated with CSS and GPS 
objects exhibit a variety of line profiles, suggesting signifi- 
cant, sometimes complex, gas motions, while the Hi column 
densities have been found to be anti-correlated with source 
sizes (Vermeulen et al. 2003; Pihlstrom et al. 2003; Gupta 
et al. 2006). Pihlstrom et al. have explored different distri- 
butions of the Hi gas density to explain this observed anti- 
correlation, and find that a disk as well as a spherical distri- 
bution are consistent with it. To study some of these aspects 
and to test for consistency with the unification scheme, we 
examine the dependence of Hi column density on relative 
prominence of the core, f c , which is being used as a statisti- 
cal measure of the orientation of the source axis to the line 
of sight (Orr & Browne 1982; Kapahi & Saikia 1982). In this 
paper, we have investigated this using a sample of 32 CSS 
and GPS sources, and have also examined the dependence of 
Hi column density and velocity shift of the primary absorp- 
tion component relative to the systemic velocity, on source 
symmetry parameters. 
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Table 1. Sample of CSS and GPS sources. 
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The columns are as follows: Cols. 1 and 2: source name and an alternative name; cols. 3 and 4: optical identification and redshift; col. 5: 
radio luminosity at 5 GHz; col. 6: the largest linear size (LLS) defined to be the separation between the outermost peaks of radio 
emission; col. 7: the fraction of emission from the core, f c , at an emitted frequency of 8 GHz; col. 8: structural classification where D 
denotes a double source while T implies a triple source which has a detected core; col. 9: the separation ratio, Rg , defined to be the 
ratio of the separation of the further component from the core to the nearer one; col. 10: the misalignment angle, A, defined to be the 
supplement of the angle formed at the core by the outer hotspots; col. 11: references for the radio structure; col. 12: radio identification 
of the sources as CSS and GPS objects; col. 13: Hi column density or a 3cr upper limit to it estimated assuming a spin temperature of 
100 K and full coverage of the background source by the absorber; col. 14: the shift of the primary Hi component relative to the 
systemic velocity as measured from the optical emission lines, with a negative sign indicating a blue-shift, and col. 15: references for the 
Hi observations. 

References for structural information: 1: Tzioumis et al. (2002); 2: Pearson & Rcadhead (1988); 3: Owsianik ct al. (1998); 4: Zensus et 
al. (2002); 5: Baum ct al. (1990); 6: Giovannini ct al. (2001); 7: Feng et al. (2005); 8: Gupta ct al. (2005); 9: Fanti ct al. (1989); 10: 
Saikia et al. (1995); 11: Sanghcra ct al. (1995); 12: Fey & Chariot (2000); 13: Dallacasa ct al. (1995); 14: Saikia ct al. (2001); 15: 
Stanghellini et al. (2001); 16: Akujor ct al. (1991); 17: Saikia ct al. (1998); 18: Taylor et al. (1994); 19: Peck & Taylor (2000); 20: 
Fomalont ct al. (2000); 21: Bondi ct al. (1998); 22: Strom ct al. (1990); 23: Liidkc ct al. (1998); 24: Mutel et al. (1981); 25: Fey et al. 
(1996); 26: Pcrlman ct al. (1996); 27: Xiang ct al. (2002); 28: Ojha ct al. (2004); 29: Xu ct al. (1995); 30: Taylor & Vermeulen (1997); 
31: Fomalont et al. (2003); 32: Polatidis ct al. (1995); 33: Taylor et al. (2000). 

References for Hi observations: V: Vermeulen et al. (2003); C: Carilli ct al. (1998); Gu: Gupta et al. (2006); PI: Peck et al. (2000); M: 
Morganti et al. (2001); P2: Peck ct al. (1999). 
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Figure 1. Neutral hydrogen column density as a function of core 
fraction. The solid symbols (in red) correspond to galaxies while 
open symbols (in blue) denote quasars. Circles denote GPS ob- 
jects while triangles denote CSS sources. Upper limits are marked 
by arrows. The vertical dashed line marks the median value of f c 
for galaxies. 



2 THE SAMPLE OF SOURCES 

Our sample consists of all CSS and GPS sources identified 
with radio galaxies and quasars for which we could compile 
information on the Hi column density and core flux den- 
sity or an upper limit to it either from the literature or 
our observations. CSS and GPS sources such as J1347+1217 
(4C+12.50) and J1407+2827 (Mrk668, OQ208), which have 
prominent Hi absorption lines and well-defined radio struc- 
tures have not been included because they are classifed as 
Seyferts in the literature (e.g. Hurt et al. 1999; Gallimore 
et al. 1999). The radio source J1415+1320 (OQ122) which 
has an Hi absorption line (Carilli et al. 1992) but has been 
identified as a BL Lac object possibly residing in a spiral 
host galaxy (Wurtz, Stocke & Yee 1996) has also not been 
considered. For this study, a core is defined to be a compact 
feature between the two outer lobes of emission. For those 
without a detected core, the upper limits to the core flux 
density are 3 times the rms noise in the image, and have 
been estimated only for those sources whose largest angular 
size is at least 5 times the resolution element. The fraction of 
emission from the core, f c , or an upper limit to it at an emit- 
ted frequency of 8 GHz has been estimated for each source. 
To estimate this fraction, we have used a spectral index of 
and 1 for the core and extended emission respectively, unless 
information was available in the literature. 

The sample consisting of 32 objects of which 14 are GPS 
and 18 are CSS objects is listed in Tabled 23 of these are 
associated with galaxies and 9 with quasars. The projected 
linear size of the sample ranges from 0.02 to 11.3 kpc with a 
median value of 0.68 kpc, while the 5-GHz luminosity ranges 
from 1.3xl0 25 to 1.2xl0 28 WHz" 1 with a median value of 
7.6 xlO 26 WHz" 1 . The values of f c range from <0.001 to 
0.279 with a median value of 0.011, while those of N(Hl) 



ranee from <0.06xl0 20 to 84.2 xlO 20 cm with a median 
value which is less than 1.8 x 10 20 cm -2 . In this paper, we will 
refer to sources with core prominence lesser and greater than 
the median value as having lower and intermediate values. 
The sources with the highest values of f c are not represented 
in our sample since they would not be classified as CSS and 
GPS objects. 



3 DEPENDENCE OF HI COLUMN DENSITY 
ON CORE PROMINENCE 

To examine the dependence of Hi column density on core 
prominence and also check for consistency with the unifi- 
cation scheme for radio galaxies and quasars, we plot the 
Hi column density against the core fraction (Fig. . In our 
sample, there is only one Hi detection for the 9 quasars and 
15 detections for the 23 galaxies, consistent with orientation 
effects playing a significant role in the detection of Hi ab- 
sorption. The quasars have a higher median value of f c of 
~0.028 compared with 0.011 for the galaxies. These trends 
are broadly consistent with the expectations for the unifica- 
tion scheme. 

One of the striking features seen in Fig. Q is that 10 
of the 15 galaxies that have a detected Hi absorption fea- 
ture, also have a detected radio core, compared with only 
2 out Of 8 sources without an Hi detection. This suggests 
that detection of a radio core in galaxies enhances the prob- 
ability of detecting a 21-cm absorption feature. Dividing the 
galaxies into two roughly equal groups at the median value 
of f c ~0.011, 9 of the 11 with f c ^0.011 (i.e. intermediate 
values of core fraction) have a detected absorption feature 
compared with 6 out of 12 detections in the other group 
with lower values of f c . The median value of N(Hl), using 
upper limits as their values, for the intermediate f c group 
is ~3.7xl0 20 cm" 2 compared with 1.5xl0 20 cm" 2 for the 
lower f c group. A Kolmogorov-Smirnov test shows the dis- 
tributions to be different at a significance level of ~95 per 
cent. Thus for galaxies, there is a tendency for the detection 
rate as well as the column density to increase with the de- 
gree of core prominence. For the 9 quasars in our sample, 
there is only one detection and hence it is not meaningful 
to do a similar exercise. But clearly, sources with the lowest 
values of f c do not tend to have the higher values of N(Hl). 

It is important to enquire whether such a trend is con- 
sistent with the unification scheme for radio galaxies and 
quasars. The deficit of high column density values for sources 
with low values of f c can be understood in a scenario where 
the radio source sizes are larger than the scale of the circum- 
nuclear disk which contributes most of the absorbing gas, so 
that the line of sight to the lobes at very large inclinations 
(corresponding to small values of f c ) does not intersect the 
Hi disk. For intermediate values of f c , one would then ex- 
pect higher detection rate as well as higher values of N(Hl) 
as seen in Fig.Q For angles close to the line of sight, which 
would have the highest values of f c , one would expect small 
values of N(Hi) and much lower detection rates due to the 
line of sight passing largely through the ionization cone and 
also through regions of the disk with a high spin temperature 
due to the presence of an AGN (Bahcall & Ekers 1969) . Such 
sources will usually get classified as compact flat spectrum 
(CFS) sources and not many of them with similar properties 
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Figure 2. Neutral hydrogen column density as a function of pro- 
jected linear size of the radio sources. See caption of Fig. for 
the meaning of the various colours and symbols. 
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Figure 4. The plot of velocities for the principal Hi absorption 
components with respect to the systemic velocity estimated from 
optical lines as a function of symmetry parameters. See caption 
of Fig. Q] for the meaning of the various colours and symbols. 
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Figure 3. Neutral hydrogen column density as a function of radio 
source symmetry paramters i.e. the separation ratio and misalign- 
ment angle. See caption of Fig. for the meaning of the various 
colours and symbols. 



as the sources in our sample have been searched for 21-cm 
absorption (cf. Gupta et al. 2006) . 



3.1 Dependence on linear size 

The neutral hydrogen column density as a function of linear 
size for our sample of 32 radio sources is shown in Fig. |5] 
This shows the anti-correlation between the Hi column den- 
sity and source size reported by Pihlstrom et al. (2003) and 
extended to a larger sample of 96 sources by Gupta et al. 
(2006). Thus detection probability of 21-cm absorption fea- 
ture is higher for both small sources and those with inter- 
mediate values of core prominence. It may be noted that 
projected linear size and f c for this sample show no signifi- 
cant relation. 



4 HI AND SYMMETRY PARAMETERS 

In addition to the circumnuclear disk, absorption could also 
occur in neutral clouds that are part of the gas that has been 
accelerated by interaction with a radio jet. Gas outflows as 
fast and broad as ~1000 km s _1 have been detected in a 
number of radio sources in neutral as well as ionised gas 
(see Morganti, Tadhunter & Oosterloo 2005). The kinemat- 
ical and spatial properties of these suggest that they are 
being driven by the same mechanism i.e. interaction with 
the expanding radio jets. The dynamical interaction of the 
radio source with these clouds could distort the morphology 
of the radio source thereby affecting the symmetry parame- 
ters such as the separation ratio and the misalignment angle. 
The separation ratio, Re, is defined to be the ratio of the 
separation of the further component from the core to the 
nearer one, while the misalignment angle, A, is defined to 
be the supplement of the angle formed at the core by the 
outer hotspots. Such jet-cloud interactions may also affect 
the velocity of the absorbing cloud relative to the systemic 
velocity of the galaxy. However, we find no evidence of any 
significant dependence of the Hi column density as well as 
the relative velocity of the primary absorption component on 
either the separation ratio or the misalignment angle (Figs. [3] 
andQ. But this does not rule out the importance of the role 
played by jet-cloud interactions in the observed structures of 
these sources as there exists significant evidence for this in 
terms of structural, brightness and polarization symmetry 
parameters (Saikia et al. 1995; 2001; Saikia & Gupta 2003). 



5 SUMMARY 

For a sample of CSS and GPS objects, we have examined 
the dependence of Hi column density on the degree of core 
prominence and projected linear size. In our sample, 15 out 
of 23 galaxies exhibit 21-cm absorption as against 1 out of 
9 quasars, which is broadly consistent with the unification 
scheme. The detection rate as well as the Hi column density 
appears to be larger for galaxies with more prominent cores, 
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with the higher column densities occurring in intermediate 
values of f c and not in the lowest values of f c . This broad 
trend can be understood in the unification scheme if the lines 
of sight towards the lobes of sources with the lowest values 
of f c , which are inclined at the largest angles to the line of 
sight, do not pass through the circumnuclear Hi disk. This 
could happen if the source sizes are larger than the scale of 
the disk. It is interesting to note that 10 of the 15 galaxies 
with a detected Hi absorption feature also has a detected 
radio core, compared with 2 out of 8 galaxies without any 
Hi detection. Our sample also shows an inverse relationship 
between the Hi column density and the projected linear size, 
suggesting that both small sizes and more prominent cores 
enhance the probability of detecting Hi absorption in CSS 
and GPS objects. 

In addition to the disk, Hi absorption may also be due to 
clouds of cold gas, some of which may have interacted with 
the radio jet. We, however, find that the Hi column density 
and the velocity shift of the primary absorption component 
relative to the systemic velocity show no dependence on the 
degree of misalignment and the separation ratio. 
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